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ABSTRACT 

We have analyzed high-resolution X-ray spectra of a sample of 22 active stars 
observed with the High Energy Transmission Grating Spectrometer on Chandra 
in order to investigate their coronal plasma density. Densities where investi- 
gated using the lines of the He-like ions O VII, Mg XI, and Si XIII. Si XIII 
lines in all stars of the sample are compatible with the low-density limit (i.e. 
rig < 10^^ cm~'^), casting some doubt on results based on lower resolution EUVE 
spectra finding densities > 10^^ cm~'^. Mg XI lines betray the presence of high 
plasma densities up to a few 10^^ cm~^ for most of the sources with higher X-ray 
luminosity (> 10'^° erg/s); stars with higher Lx and Lx/i^boi tend to have higher 
densities at high temperatures. Ratios of O VII lines yield much lower densities 
of a few 10^'^ cm~^, indicating that the "hot" and "cool" plasma resides in phys- 
ically different structures. In the cases of EV Lac, HD 223460, Canopus, /i Vel, 
TY Pyx, and IM Peg, our results represent the first spectroscopic estimates of 
coronal density. No trends in density-sensitive line ratios with stellar parame- 
ters effective temperature and surface gravity were found, indicating that plasma 
densities are remarkably similar for stars with pressure scale heights differing by 
up to 3 orders of magnitude. Our findings imply remarkably compact coronal 
structures, especially for the hotter (~ 7 MK) plasma emitting the Mg XI lines 
characterized by coronal surface filling factor, fugxi, ranging from 10~^ to 10~^, 
while we find /ovii values from a few 10"'^ up to ~ 1 for the cooler (~ 2 MK) 
plasma emitting the O VII lines. We find that /ovii approaches unity at the 
same stellar surface X-ray flux level as solar active regions, suggesting that these 
stars become completely covered by active regions. At the same surface flux level, 
/ivigxi is seen to increase more sharply with increasing surface flux. These results 
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appear to support earlier suggestions that hot 10^ K plasma in active coronae 
arises from flaring activity, and that this flaring activity increases markedly once 
the stellar surface becomes covered with active regions. Comparison of our mea- 
sured line fluxes with theoretical models suggests that significant residual model 
inaccuracies might be present, and in particular that cascade contributions to 
forbidden and intercombination lines resulting from dielectronic recombination 
might be to blame. 

Subject headings: Stars: activity — stars: coronae — stars: late-type — Sun: 
corona — X-rays: stars — plasmas 

1. INTRODUCTION 

X-ray and EUV spectra represent the most effective means to investigate the physics of 
the hot magnetically confined plasma in the outer atmospheres of late-type stars. The aim 
of spectroscopic diagnosis of these astrophysical plasmas is to probe the physical conditions 
of the gas and to derive relevant parameters such as electron density, thermal structure, 
and relative element abundances: such information should help us to constrain models of 
coronal heating and structuring. As an intermediate step toward this goal, by studying 
the emission properties of stellar coronae we can test the validity of the solar analogy, i.e. 
whether the basic components of stellar coronae are analogous to those observed in the solar 
corona, whose dominant X-ray emitting structures are defined by the loop-like morphology 
of coronal magnetic fields. 

The new generation X-ray observatories, Chandra and XMM-Newton, with their un- 
precedented spectral resolution and large effective providing us with high-resolution 
coronal X-ray spectra in which individual spectral lines can be readily resolved. The X-ray 
plasma diagnostics these spectra now allow for stars were previously accessible only for the 
Sun. It is also worth noting that Chandra and XMM-Newton stellar coronal spectra offer 
significantly more extensive X-ray wavelength coverage than achieved with solar instrumen- 
tation. In particular, with regard to the paper in hand, the High-Energy Transmission 
Grating Spectrometer (HETGS) on C/iandra provides high resolution (A/AA ~ 100 — 1000) 
spectra in the energy range 0.4-8.0 keV (1.5-30A) where we observe simultaneously a large 
number of prominent emission lines that, at least in principle, provide both temperature and 
density diagnostics. 

Of special importance for the purposes of this paper are the complexes of He-like 
"triplets" of O, Ne, Mg, and Si, which include the resonance line (r : Is^ ^Sq - ls2p ^Pi), 
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the intercombination lines {i : Is^ ^Sq - ls2p ^P2,i), and the forbidden hne (/: Is^ "^Sq - ls2s 
^Si), which correspond to transitions between the n = 2 shell and the ground level, n = 1. 
The utility of these lines for plasma diagnostics was pointed out by Gabriel & Jordan (1969) 
(see also Pradhan & ShuU 1981, and Porquet & Dubau 2000 for a recent treatment): the 
ratio R = f /i is mainly density sensitive, owing to the metastable ls2s ^Si level, while the 
ratio G = {f + i)/r is mostly temperature sensitive. Similar diagnostic techniques have been 
widely exploited for electron density estimates in EUV spectra (see, e.g., Mariska 1992 for so- 
lar applications, and. Laming 1998 for a review of stellar work). Chandra and XMM-Newton 
data offer the opportunity to apply these techniques to the X-ray range and to explore dif- 
ferent temperature and density regimes. In this context, analyses of He-like triplets seen 
in the X-ray spectra of a few different stars have appeared in the recent literature (e.g., 
Brinkman et al. 2000; Canizares et al. 2000; Audard et al. 2001; Giidel et al. 2001ab,2002; 
Huenemoerder et al. 2001; Ness et al. 2001,2002; Raassen et al. 2002; Stelzer et al. 2002, 
Sanz-Forcada et al. 2003b). The study of Ness et al. (2002) is more extensive and involves 
several stars observed with the Chandra Low Energy Transmission Grating Spectrograph 
(LETGS). 

While all the results from these studies cannot be easily summarised, it is apparent, 
mostly from O VII lines, that there is a general trend showing that plasma with temperatures 
of up to a few million K has densities similar to that found in solar active regions — up to a 
few 10^° cm~^ or so. Based on lines of Mg XI, there is also some evidence for higher densities 
around ~ 10^^ cm~^. These results, however, are not without controversy. Brickhouse (2002), 
for example, notes that different density values have resulted from different analyses of the 
same stars observed with the different Chandra and XMM-Newton spectrometers. These 
instruments have different resolving powers, and so density-sensitive lines can be blended to 
greater or lesser extents with other lines. In this regime, the analysis is also sensitive to the 
continuum level adopted. For the case of XMM-Newton spectra, this is further complicated 
by the extended wings of the instrumental profile, which can give rise to a significant pseudo- 
continuum. 

The best instrumentation to study the He-like lines of O, Ne, Mg and Si is that with 
the highest resolving power — the Chandra HETGS. We were therefore motivated to perform 
a systematic study of these lines in a large sample of stars observed with the HETGS in 
order to attempt to characterize plasma physical conditions in a sample of stars with very 
different stellar parameters. The results can be used to probe possible trends of the plasma 
parameters with more global coronal parameters, such as X-ray luminosity or flux, or stellar 
effective temperature, gravity, rotational period, etc. 

In this paper, we have analyzed in detail the lines of Si XIII (6.7A), Mg XI (9.2A), 
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and O VII (22A), in a sample of 22 late-type stars, covering a range of spectral types and 
activity levels. We did not undertake an analysis of the Ne IX triplet lines here, since they 
are heavily blended with iron (predominantly Fe XIX; see Ness et al. 2003). The Si, Mg and 
O He-like triplets offer a reasonable sampling of different temperature and density regimes, 
with temperatures of formation and density sensitivities approximately corresponding to the 
following values: 10 MK, lO^^ cm'^ (Si XIII), 7 MK, 10^^ ^j^-s (^^g XI); 2 MK, 10^° cm'^ 
(O VII). 

We describe the observations and the stellar sample in §2, and the detailed line analysis 
in §3; we discuss the results obtained in the context of coronal activity in section §4 and we 
draw our conclusion in 55. 



2. OBSERVATIONS 

The sample analyzed consists of all the late-type stars for which observations were 
available in the Chandra public archive^ at the time when this study was undertaken. The 
sample comprises 22 cool stars of different activity levels, and includes single stars (the flare 
stars AU Mic, Proxima Cen and EV Lac, the zero-age main-sequence star AB Dor, the 
T Tauri star TW Hya, and the giants HD 223460, 31 Com, (3 Get, Canopus and fi Vel) and 
multiple systems (the close eclipsing binary system Algol; two eclipsing binary systems of 
W UMa type, namely ER Vul, and 44 Boo; and the RS CVn systems TZ CrB, UX Ari, 
^ UMa, II Peg, A And, TY Pyx, AR Lac, HR 1099, IM Peg). 

These sources were all observed using the Chandra-HETGS and the Advanced CCD 
Imaging Spectrometer (ACIS-S) detector (Canizares et al. 2000). The High-Energy Trans- 
mission Grating disperses X-rays through two different gratings, whose spectral coverage 
partially overlap; they are optimized for "medium" (MEG) and "high" (HEG) energies. The 
spectral resolution obtained with the HETGS is in the range A/AA = 100 — 1000. The 
data used here were either recently processed by the standard pipehne and obtained directly 
from the archive, or else, in the case of older observations, reprocessed using standard CIAO 
3.0 tools and analysis threads to take advantage of recent improvements in the spectrum 
reduction software and ancillary calibration data. Our analysis of spectral lines was based 
on the summed data from positive and negative orders. 

Table 1 and Table 2 summarize, respectively, the stellar parameters and the charac- 
teristics of the observations. For each source the listed X-ray luminosity was derived by 
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integrating the photon energies over all the wavelengths in the spectral ranges of the two 
different instruments, HEG and MEG. 

As summarized in Table 1, the single stars in our sample are four late-type flare stars, 
namely AU Mic, Proxima Centauri, EV Lac, and the rapid rotator AB Dor (Prot ~ 44 ksec), 
the giants, HD 223460, 31 Com, j3 Get, Ganopus and fi Vel and the T Tauri star TW Hya. 
The other sources are binaries and most are classified as systems of RS GVn-type. In detail, 
among the RS GVn systems: II Peg is a system of a K2 V star and an unseen companion; 
analogously A And is a binary system of which we see only one component (G8 IV-III); 
finally, TY Pyx is an eclipsing active close binary system of two G5 IV stars. As for the 
remaining sources, ^ UMa is a binary system of two optically identical components (GO V) 
one of which is a rapid rotator, and there are two eclipsing binaries of W UMa type: ER Vul, 
a short period active binary system of solar-like stars (GO V + G5 V) that are rapid rotators 
(v~ 81, 71 km/s respectively), and 44 Boo, a contact binary system of G1V-G2V stars. 

It is worth noting that our sample contains sources whose X-ray luminosities span a 
very wide range — from a few 10^^ erg/sec of Proxima Gen up to a few 10^^ erg/sec of the 
giant HD 223460. As for the ratio of X-ray to total luminosities, Lx/i^bob our sample covers 
the range Lx/Li^oi ~ 10~^-10~^. We are therefore analyzing spectra of stars whose activity 
levels cover almost the whole activity range observed for late-type stars, although most of 
them are very active stars. 

Spectra for each source were extracted from photon event files obtained from the Ghan- 
dra Data Archive^ using GIAO 3.0. Effective areas were calculated using standard GIAO 
procedures, and we included the effects of the AGIS-S contamination layer (e.g. Plucinsky 
et al. 2002) using Version 1 of the effective area contamination correction^. The result- 
ing HETG spectra are shown in Figures 1 and 2; spectra of sources of the same type are 
grouped together: Fig. 1 shows the spectra of single stars (both dwarfs and giants); Fig. 2 
shows the spectra of binary systems. For sources with several different observations we show 
the coadded spectra. 
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3. DATA ANALYSIS 
3.1. Lightcurves 

As a first step of the analysis we have extracted the lightcurves to check whether the 
spectra integrated over the whole observation time are influenced by significant variability 
or large flare events. In order to obtain the lightcurves we excluded the 0th order events that 
can be drastically compromised by photon pile-up effects, and considered just the dispersed 
events. The lightcurves were obtained by summing all the counts of HEG and MEG spectra 
for spectral orders n > 1 in temporal bins of 100s. 

The resulting lightcurves are shown in Figures 3 and 4. We see from these figures that 
there are no remarkable variations in X-ray count rate for 31 Com, (3 Ceti, Canopus, /i Vel, 
ER Vul, UX Ari, ^ UMa, A And, AR Lac, and IM Peg. In the remaining cases, the lightcurves 
of AU Mic and EV Lac show flaring activity, and in particular, flares of different intensity 
and temporal extent; the Proxima Centauri lightcurve shows a flare at the beginning of the 
observation and then a decay over about 3 ksec and then the emission stabilizes at a low 
level; the AB Dor emission is almost steady with the exception of a moderate flare extending 
over about 10 ksec (of the ~ 53 ksec of the total observation); HD 223460 shows a large flare 
over half of the observing time — it is worth noting that flares have been thought of as fairly 
unusual on single giants (e.g. Ayres et al., 1999,2001); the lightcurve of Algol shows evidence 
of a big flare with a peak at about 10 ksec after the beginning of the observation and then a 
slow decrease over about 40 ksec to approximately the initial emission level; TZ CrB shows 
a flare with a rapid increase of the emission by a factor 4 over the last ~ 3.5 ksec among the 
83.7 ksec of total exposure time; II Peg shows a big flare analyzed in detail by Huenemoerder 
et al. (2001); the TY Pyx total emission is quiescent with a slight modulation for almost the 
whole observation except a short flare (lasting ~ 3 ksec) in the middle of the observation. The 
lightcurve of 44 Boo exhibits periodic modulation that was analyzed in detail by Brickhouse 
et al. (2001) who exploited, for the first time, the possibility of a crude Doppler imaging at 
X-ray wavelengths to derive geometrical properties of coronal X-ray emitting structures. 

Regarding the influence of these variations on the integrated spectra, in most of the cases 
they do not affect substantially the basic properties of line emission, whereas the continuum 
emission changes significantly during the flares. For example, for the II Peg observation, 
Huenemoerder et al. (2001) show that at "low" temperatures (T< 10 MK) the modulation 
of the line fluxes during the flare is low (< 15%), and in particular, since the triplet lines do 
not show substantial variations, they use the data integrated over the entire observation for 
the density diagnostics. In other cases, e.g., AU Mic, TZ CrB, TY Pyx, AB Dor, the flares 
extend over a very small portion of the observing time, thus they do not affect substantially 
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the total spectra. Unfortunately, in these cases the flaring intervals are also too short to 
provide adequate signal for separate analysis. Finally, in the case of Proxima Centauri, the 
flaring portion of the emission is the only detectable signal, and even when considering the 
total spectra, it is not possible to measure the fluxes of most of the lines because of the 
low signal. In the light of the above considerations, our analysis here is based on spectra 
integrated over entire observation. 

3.2. Spectra 

We illustrate the lines used in this analysis in the MEG spectrum of (3 Ceti, illustrated 
in Figure 5; the relevant lines are labeled by element and ionization stage. As already men- 
tioned in Section 1, in the HETGS spectral range the lines of He-like ions are among the 
most intense, though our analysis hinges on the measurement of relatively weak components 
of these. The theoretical wavelengths of these line complexes, as they appear in the APED 
database (Smith et al. 2001), are listed in Table 3. Our plasma density analysis basically con- 
sists of a comparison between measured and theoretical relative intensities of the forbidden 
and intercombination lines (the R ratio). For this comparison, we adopt the line emissivities 
of APED, and in Figure 6 we show the dependence of the R ratio on at the temperature 
of maximum formation of each of the He-like ions considered here. For comparison we plot 
also the values from Porquet et al. (2001). 

As Figure 6 illustrates, the R ratios at the low density limit differ by about 10% between 
the different calculations. At densities above the lower sensitivity threshold, the calculations 
are in much better agreement for O and Mg, though differences persist for Si at higher 
densities. The Porquet et al. (2001) and Smith et al. (2001) emissivities are based on different 
input atomic data, and so small discrepancies are to be expected. Mewe et al. (2003) noted 
the same differences in predicted R ratios for O VH, and attributed these to the different 
sources used for the direct excitation terms: Smith et al. (2001) used the excitation rate 
coefficients from Kato & Nakazaki (1989) (for n< 5) and from HULLAC (Liedhal et al. 
1995; for n=6-10), while Porquet et al. (2001) adopted coefficients from Zhang & Sampson 
(1987) (for n< 3) and from Sampson et al. (1983) (for 2 < n < 6). We also note that 
dielectronic recombination and radiative recombination rates are from different sources in 
the two models. Since the 1525^5"! and \s2p^Pi 2 levels are primarily populated by cascades 
from higher levels that in turn are largely populated by dielectronic recombination (e.g. Smith 
et al. 2001), it seems plausible that the treatment of these processes could also contribute 
significantly to differences in predicted R ratios. 

In this paper, we emphasise as much as possible trends in the observed R ratios, rather 
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than in the exact conversion to density values using theoretical models; had we adopted the 
Porquet et al. (2001) theoretical R ratios instead of those of Smith et al. (2001), none of our 
conclusions would be significantly changed. 

Observed spectra were analyzed with the IDL 5.3^ PINTofALE 1.0^ software (Kashyap 
& Drake 2000). Line fluxes were determined by fitting with a modified Lorentzian function 
described by the relation 

F{X)=a/il + i^)r (1) 

where a is the amplitude and F is a characteristic line width, and with a value of 2.5 for the 
exponent j3. This function has been found to be a good match to observed line profiles in 
stellar coronal spectra with counts of a few thousand or less (see, e.g., the Chandra Proposers' 
Observatory Guide^). 

The measured R ratio can be somewhat sensitive to the level assumed for the underlying 
continuum, and so we took particular care in locating this. For a few spectra of the very active 
stars the continuum is quite strong in the vicinity of the Mg and Si triplets (HD 223460, Algol, 
TZ CrB, II Peg, AR Lac, HR 1099, and IM Peg), while in the 21-22 A range where the O VII 
lines are located the continuum level was always close to zero and essentially negligible. The 
continuum in a coronal plasma at temperatures exceeding a few million K is predominantly 
due to bound-free and free-free emission from hydrogen. Over small wavelength intervals 
this continuum is quite fiat and can be modelled using the instrument effective area function 
pinned down in adjacent "line-free" regions. For the very active stars, we also used a generic 
continuum model computed for the RS CVn-type binary HR 1099 by Drake et al. (2001), 
scaled again according to spectral regions deemed essentially "line-free" . We also performed 
sensitivity tests to determine how errors in continuum location effected our measured fluxes: 
line fluxes were measured for continuum levels deemed too high and too low, based on our 
location criteria described above. In the case of the / and i lines, the differences in measured 
fluxes were always smaller than the statistical uncertainties in the fluxes themselves so that 
uncertainty in continuum placement was not a significant source of error in the analysis. 

We measured the spectral line intensities of the Si XIII, Mg XI and O VII He-like 
triplets from both HEG and MEG spectra wherever the data offer useful signal. The results 
obtained from the two spectra are in good agreement, as we will show below. The O VII 
lines were measured only in the MEG spectra since they are outside the wavelength range 
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of the HEG. For all the He-like line complexes we fitted the three lines simultaneously; the 
free parameters of the fit were the position and the flux of each of the lines, while the widths 
of all lines were constrained to have the same value. 



3.3. Mg XI spectral region 

While the Si XIII and the O VII triplet lines are relatively isolated, the Mg XI triplet lies 
in a more crowded region and the lines are more likely to be affected by blending. Indeed, the 
APED database predicts a number of Fe XIX-XXII lines in the 9.15-9.35 A region, though 
each of these is considerably weaker than the Mg lines of interest. Moreover, lines from the 
Lyman series of H-like Ne with upper levels n > 5 also lie in this region. These latter lines 
are listed in Table 4. 

We investigated the possible influence of line blends on our measurement of the Mg XI 
triplet components in some detail. First, we carefully examined HEG spectra with both 
the highest signal-to-noise ratio and the highest relative intensity of Fe or Ne lines with 
respect to Mg lines (e.g., UX Ari for both Fe and Ne, AB Dor for Fe, HR 1099, II Peg and 
AU Mic for Ne) for evidence of spectral features close to, or blending with, the triplet lines. 
In particular, in the spectra of stars in which Ne has been found to be strongly enhanced 
relative to Fe, such as HR 1099 (Drake et al. 2001) and II Peg (Huenemoerder et al. 2001), 
the Ne Lyman series lines at 9.215 (9 ^ 1), 9.246 (8 1) and 9.291 (7 1) are readily 
discernible. Some lines of highly ionised Fe also probably blend with the Ne features at 
some level, though the dominant components appear fairly well-centred on the expected Ne 
wavelengths. A simple fit to the Mg components of the HR 1099 in the HEG spectrum, 
ignoring the Ne blends, is illustrated in Figure 7 {top panel). It is clear from this figure that 
the fit to the intercombination line also encompasses a significant fraction of the Ne lines: 
ignoring these blends results in a spuriously large intercombination line flux. 

In order to account for the blends, we constructed an empirical model of the region that 
included the major Mg XI lines, together with blending components. These components 
comprised four additional lines corresponding to the positions of the Ne Lyman series lines 
with upper levels n = 7, 8, 9, 10. While it is desirable to include explicitly any significant 
lines of Fe, this presents difficulties because the wavelengths of the relevant Fe lines in this 
region are poorly known. To our knowledge, the best wavelengths currently available for this 
particular region are those predicted by the APED database, and these are based primarily 
on theoretical calculations, rather than on experiment. The model itself was constrained by 
comparison with, and by fitting to, co-added HEG spectra with the highest signal-to-noise 
ratio (/5 Ceti, TZ CrB, AR Lac and HR 1099). Since Fe features might distort and displace 
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the centroids of the Ne features, in the fitting process we allowed the wavelengths of the 
four blends to vary. Line widths were all constrained to the same value, and the wavelength 
separations of the Mg lines were also fixed at their accurately known theoretical values. 

In Figure 7 we show the best fit model together with the individual components super- 
imposed on the coadded HEG spectra. The wavelength intervals between the blends and 
the Mg XI resonance line found from this fit were then fixed into the model, which was 
then tested on the coadded MEG spectra of the same four sources. This latter fit, which is 
also illustrated in Figure 7, provided an excellent match to the MEG spectrum, and yielded 
relative line fluxes in agreement with the HEG values to within statistical uncertainties. The 
individual wavelengths of the different line and blend components in the final model of the 
Mg XI spectral region determined in this way from the coadded HEG spectra are listed in 
Table 5, together with the identifications we attribute to the features. 

While the blending components include the contribution of all the Ne lines of the Lyman 
series present in this spectral range, there are doubtless contributions at some level from Fe 
lines, and we list candidates in Table 5. The relative strengths of the Ne Lyman lines for the 
higher n transitions should scale in proportion to the oscillator strengths of the transitions. 
This appears to hold for those with upper levels n = 7,8 and 9; however, the 10 ^ 1 
transition at 9.194 A appears too strong and it is shifted slightly blueward of its predicted 
location. Examination of the APED database reveals an Fe XXI transition in this vicinity 
that we attribute to this. In comparison to this Fe XXI line, other potential Fe blends 
predicted in APED appear significantly weaker by factors of at least 3 and are likely less 
significant. The wavelengths fitted for the other blend components also remained very close 
to the values of the Ne lines, and the presence of further significant blends that we have 
omitted seem unlikely. 

We compare in Figure 8 the R and G ratios obtained for the Mg XI triplet in each star 
for the case in which only the continuum and Mg XI lines themselves were included in the 
line fitting with the ratios obtained when we fitted the region using the model including 
blends. As expected, because the intensity of i is generally diminished compared to that of / 
in the latter approach, the effect is of raising the R ratio and slightly decreasing the G ratios 
compared to the former approach, i.e. including blends implies shifts toward lower derived 
densities and slightly higher temperatures. 

As an additional test that our final Mg line fluxes, and especially those of the density- 
sensitive i and / pair, are not significantly affected by unseen blends, we have investigated 
the observed f /i ratios as a function of both Ne and Fe relative line strength. As a measure 
of relative line strengths of Ne and Fe, we chose the ratios of the Mg XI resonance line with 
Fe XXI 12.29 A, Fe XVIII 14.21 A, and Ne X 12.13 A. The two lines of iron at different 
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ionization stages provide tests at different temperatures, which could be important since the 
thermal structure of our sample stars that cover a fairly wide range in X-ray activity might 
be quite varied. The observed R ratios are illustrated as a function of our control ratios in 
Figure 9. 

We note from Figure 9 that the relative intensities of the chosen Fe and Ne lines with 
respect to the Mg resonance line span a wide range of factors of 4-5 (Fe XVIII,XXI) and 
~ 15 (Ne X). Both the thermal structure and the elemental abundances of the emitting 
coronae contribute to produce these differences; the Mg XI resonance line, for example, is 
formed at temperatures intermediate between those of Fe XVIII and Fe XXI. We might also 
expect some degree of correlation between coronal density and both temperature structure 
and chemical composition. 

Of the three panels in Figure 9, only the Mg XI R ratio variation as a function of Mg 
r/Fe XVIII line strength appears completely flat, indicating an R ratio that is independent 
of Fe XVIII line strength. This is encouraging because Fe XVIII can be considered a rough 
proxy for the coronal Fe abundance, and in the case where remaining Fe blends were affecting 
our derived R ratios, we would expect to see some correlation between the two. The hidden 
Fe lines that we might expect to blend with the Mg lines come from slightly higher ionisation 
stages, however. In the case of Fe XXI, the R ratio appears somewhat higher for the very 
lowest Fe XXI/Mg r values, with a step to lower R ratio for Fe XXI/Mg r > 1.3. While the 
reason for this is not clear — as we emphasise above, both temperature structure and relative 
abundances play a role in the Fe/Mg and Ne/Mg ratios — we discount blends as the culprit 
because the trend in R ratio remains flat while Fe XXI/Mg r varies by a factor of 3-4. The 
Fe XXI 12.29 A line is also sensitive to density, and decreases in intensity with densities 
above several 10^^ cm~^. While this could contribute to the trend in Figure 9, we do not 
expect this to be a major effect since densities of order 10^^ cm"'^ or more would be required; 
in the analysis described below we rule out such high densities based on our observed Si XIII 
R ratios. 

Perhaps the most interesting trend is that between R and Ne X/Mg r: there does appear 
to be a weak correlation here, with Ne X/Mg r increasing as the R ratio decreases. At face 
value this might suggest residual Ne blending problems; however, the Ne Lyman series lines 
have very precise wavelengths and our blending model should accurately account for them. 
We conclude that the observed trend is real: that there is probably a weak trend of increasing 
Ne X line strength with increasing R ratio. We return to this in the discussion (§4). 



- 12 - 



3.4. Results 

Table 6 and 7 summarize the results for the measured line fluxes and la uncertainties 
for all the stars. The uncertainties were obtained from a full search of parameter space 
for each spectral fit. HEG and MEG line fluxes were found to differ beyond their 1 a 
uncertainties to a small degree. We first investigated to determine whether or not this could 
be the result of the spectral resolution difference between them, resulting in systematically 
different estimated continua, but we were able to rule this out: the required continuum 
placement errors are implausibly large. Instead, we have found that similar flux differences 
are obtained by integrating spectral flux in arbitrary intervals across the bandpass. We 
attribute this to effective area calibration errors, which seem to be at the 10-20% level at 
this time.^ Flux differences across the bandpass vary only very slowly with wavelength such 
that over narrow intervals the effect is grey and does not influence the measured line ratios. 

Figure 10 shows examples of the HEG spectra in the region of the Si XIII He-like triplet 
for one star of each class: single stars, giants and multiple systems. We plot the different 
components of the best fltting model (continuum emission and single spectral line features), 
together with their sum, superimposed on the data. 

Analogously, Figure 11 and Figure 12 show examples of the MEG spectra for the Mg XI 
and O VII He-like triplet regions, together with spectral flts. None of the spectra obtained 
for giant stars yielded signiflcant O VII line measurements. The spectrum of TW Hya is 
extremely interesting, since the forbidden line is essentially not detected, and thus it is the 
only case of line fluxes at the high density limit for oxygen (see also §4 and Kastner et al. 
2002). 

Figure 13 summarizes all the results obtained for the R ratios for all the sources, plotted 
as a function of surface X-ray flux and X-ray luminosity. Some of the ratios, mostly from 
HEG spectra, are missing due to the low S/N ratio of the relevant features (see Table 6). 
We use different symbols for different classes of sources as indicated in the flgures. It is 
clear from these flgures that HEG and MEG measurements are generally compatible within 
experimental uncertainties. The dashed lines of Figure 13 mark the low-density limits of 
the sensitivity of the different diagnostics, i.e. R ratios consistent with the dashed lines to 
within experimental uncertainty provide only an upper limit to the plasma density, while 
those signiflcantly below the limit provide density estimates (see Fig. 6). 

At face value, all the Si XIII R ratios are compatible with the low-density ratio limit. 
He ^ 10^^ cm^. However, we note that a good fraction of the measured ratios lie above this 
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limit. This suggests that either the theoretical low density limit R ratio is too low, or that 
some fraction of our observed R ratios are spuriously high. We return to this in §4. 

For the Mg XI triplet lines, formed at lower temperatures with respect to the Si XIII 
lines, in many cases the R ratio is lower than the low-density limit value and yields useful 
density diagnostics. The Mg XI R ratio vs. surface flux (Figure 13) shows a discernible 
trend that stars yielding high densities tend to have higher fluxes. One exception is the 
low surface flux supergiant Canopus (Fx = 0.17 x 10^ erg/cm^/s), whose Mg triplet lines 
point to high density. An analogous trend appears in the figures illustrating the R ratios 
vs. X-ray luminosity; in particular for the Mg XI R ratio this trend with luminosity is even 
more pronounced than with surface flux. One exception is the late-type single star EV Lac, 
which, owing to its late spectral type, has a relatively low X-ray luminosity but apparently 
high plasma density. 

Since the He-like triplet features are mildly sensitive to temperature, in order to derive 
actual density values from the measured R ratios we need to determine the plasma tempera- 
tures at which the density- sensitive lines are formed. The natural choice would be to use the 
temperature-sensitive G ratio; however, based on our measurements the G ratio seemed to 
yield temperatures systematically lower than what we would expect. The true temperature 
of line formation depends on both the line emissivity as a function of temperature (that 
peaks at ~ 6.8MK in the case of Mg XI and ~ 2MK for O VII), and on the thermal struc- 
ture of the corona, i.e. the emission measure distribution; since most of these active stars 
have emission measure distributions peaked around 10 K or even higher (e.g., Ayres et al. 
1998 for some giants, Drake et al. 2001 for HR 1099, Sanz-Forcada et al. 2002 for UX Ari 
and AB Dor, Huenemoerder et al. 2001 for II Peg, Singh et al. 1996b for AR Lac), we then 
expect to get temperatures even higher than the nominal maximum formation temperature. 
Instead, typical derived temperatures from the G ratio were lower than those of the peak 
emissivities for both Mg XI and O VII: for Mg XI, typical derived temperatures were in the 
range 1-6 x 10^ K, while for O VII temperatures where only 5 x 10^-2 x 10^ K. 

Since the G ratio temperatures seem too low, we compared them to ionization tem- 
peratures derived by comparing predicted and observed ratios of the He-like to H-like Lya 
resonance lines of O and Ne: the representative ionization temperature is that at which an 
isothermal plasma most closely matches the observed ratio. Line fluxes for the O and Ne 
Lya lines were again measured using the PINTofALE FITLINES routine, though we only 
examined the lines in MEG spectra; fluxes are listed in Table 7. The comparison between G 
ratio and ionization temperatures is illustrated in Figure 14. The G ratio yields temperatures 
systematically lower than ionization temperatures by factors of ~ 2-4. 

The scatter in the derived ionization temperatures for the different stars is very small 
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and amounts to less than a factor of 2. The ionization temperature probably overestimates 
very slightly the true mean temperature of formation of the He-like lines, since it represents 
the temperature intermediate between the H-like and He-like ions. Since we discount the G 
ratio temperatures as spuriously low, and the ionization temperatures as slightly high, for 
simplicity in the derivation of densities we adopted the temperature of peak emissivity for 
both O VII and Mg XI. It is worth emphasising here that the R ratio is only weakly dependent 
on temperature, and electron densities derived using all the different temperatures obtained 
with different methods yield density values that differ by a factor at most 2 — considerable 
lower in most cases than the statistical uncertainty associated with measurement errors in 
the R ratio. The final density values as derived from the R ratios are illustrated as a function 
of both X-ray luminosity and surface flux in Figure 15 and listed in Table 8. 



4. Discussion 

The survey described here provides information on plasma densities in the X-ray emit- 
ting coronae of a wide sample of stars. We feel in particular that our measured fluxes are 
robust, and that the deblending procedure discussed in §3.3 provides for much more accurate 
measurement of the density- sensitive Mg XI lines than has been presented in other studies 
to date. 



4.1. Atomic Data 

The derived plasma densities we will discuss below depend on the predictions of atomic 
models, and it is natural to ask what systematic errors might be associated with converting 
the measured R ratios into electron densities. In particular, our study highlights two potential 
problems with recent calculations of He-like emissivities for Si, Mg and O. 

In the case of Si, the APED predicted low density limit value of the R ratio is about 
50% lower than what we observe, assuming the source plasmas are indeed at the low density 
limit. The Porquet et al. (2001) predicted ratio is slightly higher (Figure 6), though still 
discrepant with regard to our measured ratio. Instead, for Mg the low density limit R ratio 
from both APED and Porquet et al. (2001) agrees very well with the upper envelope of the 
observed ratios, as we would expect if some of the coronae were at the low density limit 
and some comprised higher density plasma. The case of O VII is more difficult to evaluate, 
since it appears that all the measurements imply densities above the theoretical low density 
limits, though some of these are of course not statistically significant. 
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The temperatures derived from the APED G ratios of both Mg and O are much lower 
than we would expect for these coronal plasmas. In the case of O VII, the observed G ratios 
should have reached a limiting value of ~ 0.6 for most of the stars in our study, according 
to the APED emissivity ratio. The situation is improved for the case of the Porquet et al. 
(2001) ratios, though G ratio temperatures still seem too low. 

While direct excitation dominates other mechanisms in populating the ls2p ^Pi upper 
level of the resonance line, in the case of the forbidden line cascades from higher n levels are 
much more important. These cascade contributions predominantly come from dielectronic 
recombination. The G ratio is therefore quite sensitive to the balance between dielectronic 
recombination and direct excitation processes. Since our G ratio temperatures are too low, 
this suggests that either the model populations of the / and i upper levels are underestimated, 
or the resonance line, r, ^Pi upper level population is overestimated. Considering that the 
predicted Si R ratio low density limit seems too low to explain the observations, our results 
suggest that the culprit is likely in the cascade contributions and the equipartition between 
cascade rates into ^P and ^S. The size of the discrepancy we see suggests errors in the 
relative level populations of about 20% or so. As we noted earlier in connection with the 
comparison of Smith et al. (2001) and Porquet et al. (2001) predicted R ratios, uncertainties 
of at least 10% or so are clearly present in these calculations. 

In this study, we emphasize trends with the observed R ratio itself, in addition to 
the values of the derived plasma densities. While it appears that only the Si R ratio is 
significantly affected by the model discrepancies described above, it should be kept in mind 
in the following discussion that there is an inherent systematic uncertainty in converting the 
observed ratios to densities, in addition to photon counting statistics. 

4.2. Densities 

We characterise the general R ratio and plasma density findings here as follows: 

1. R ratios and, consequently, plasma densities at temperatures of ~ 10^ K are all very 
similar, with typical densities of rie ~ 2 x 10^° cm~^ and a scatter between different 
stars of approximately a factor of 2. 

2. We find a trend of higher plasma density in sources with higher coronal temperatures, 
and with higher coronal X-ray luminosities. This trend is less obvious when cast in 
terms of surface flux, and more obvious when viewed as a function of the ratio of X-ray 
to bolometric luminosities. 
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3. In all the stars studied, the observed Si XIII R ratios are above or similar to the 
predicted low density limit. At face value this indicates that no star has a coronal 
plasma density > 10^^ cm~'^ at temperatures of ~ 10^ K or higher. 

The general finding of higher plasma densities at higher temperatures in active stars 
has already been suggested by the growing bulk of work from other EUV and X-ray obser- 
vations of active stars (e.g., Bowyer, Drake & Vennes 2000 and Drake 2001 for summaries of 
EUVE studies; Brickhouse 2002, Sanz-Forcada et al. 2002,2003ab, Huenemoerder et al. 2001, 
Argiroffi et al. 2003). Our survey here confirms the findings of the earlier work; the EUV 
studies were often not unambiguous, owing to the possible influence of lines blending with 
the density-sensitive diagnostics (e.g. Drake 2001). Indeed, the upper limit of rig < 10^^ cm""^ 
imposed by Si XIII for all stars of our sample casts doubt on findings from some Fe EUV 
line ratios of densities Uc ~ 10^^ cm~^, as reported, e.g., by Sanz-Forcada et al. (2003a), 
though we have to keep in mind possible inaccuracies in the theoretical Si XIII R ratio. 

The results for different groups of stars are summarized below. Among our results, 
coronal density estimates are presented here for the first time for several sources of our 
sample: the single dwarf EV Lac, the giants HD 223460, Canopus and fi Vel, and the 
RS CVns TY Pyx and IM Peg. For the remaining sources, the results we obtain here are 
compared with findings of previous works in Tables 9 and 10. 

4-2.1. Single Dwarfs 

In single dwarfs, excluding the peculiar case of TW Hya that we will shortly discuss 
below, high density seems to be associated to the presence of evident flaring events as 
observed for AU Mic (Monsignori Fossi et al. 1996; Magee et al. 2003), Proxima Gen (Giidel 
et al. 2002), and EV Lac^ (this work). 

The results obtained for single dwarfs are generally compatible with the results present 
in the literature (see Table 9), at lower temperatures (2-3 MK), while some disagreement 
is present for higher temperatures (6-10 MK). For example for the rapidly rotating zero 
age main sequence KO dwarf AB Dor, our density estimates are in good agreement with 
the finding of Sanz-Forcada, Maggio & Micela (2003b), based on both XMM-Newton and 
Chandra spectra, though their estimate from Mg XI is slightly higher than our estimate. 



^EV Lac X-ray observations prior to the Chandra and XMM-Newton era comprise only low resolution 
studies, including those of BeppoSAX and ROSAT (e.g., Sciortino et al. 1999), and ASCA (e.g., Favata et 
al. 2000), in which direct density diagnostics are unavailable. 
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possibly because of the blending due to Ne and Fe that we have taken into account here. 
However, our results for Mg are somewhat in conflict with the analysis of Fe XIX-XXII lines 
in EUVE, Chandra and XMM-Newton spectra by Sanz-Forcada et al. (2002,2003a), who find 
somewhat higher densities of rig ~ 10^^ — 10^^ cm~^. In particular, the densities estimated 
from X-ray Fe lines by Sanz-Forcada et al. (2003a) span an order of magnitude. While it is 
possible that the slightly higher temperatures probed by these Fe lines could be characterized 
by significantly higher densities than the plasma responsible for the Mg XI lines, the limits 
placed on the density by the Mg lines tend to rule out the presence of significant amounts 
of plasma at densities close to 10^'^ cm~^. 

TW Hya The same Chandra-HETG spectra of the pre-main sequence T Tauri star TW Hya 
have been analyzed by Kastner et al. (2002), who reported similar results to those obtained 
here: the Mg XI lines have not been detected, and the O VII lines yield ric > 10^^ cm~'^. 
The O VII spectral region of TW Hya is strikingly different from the spectra of all the other 
sources in our sample that yield ric invariably close to 10^° cm~'^. Stelzer & Schmitt (2004) 
analyzed XMM-Newton observations of TW Hya at both high (with the Reflection Grating 
Spectrometer) and intermediate (with the European Photon Imaging Camera) spectral res- 
olution, finding densities based on O VII in good agreement with our results and those of 
Kastner et al. (2002). Since, during the XMM-Newton observations no clear flare is detected, 
Stelzer & Schmitt (2004) conclude that the unusually high density found in TW Hya is not 
due to any flaring activity. Echoing Kastner et al. (2002) and Stelzer & Schmitt (2004), the 
X-ray emission from TW Hya seems radically different from that typical of stellar coronae, 
supporting conjectures that it originates from a fundamentally different plasma, such as one 
heated by a shock at the bottom of an accretion column. 

4-2.2. Single Giants 

None of the spectra of the single giants analysed here yielded useful measurements of 
the O VII lines. In the case of Mg XI, however, we find clear indication of high densities of a 
few 10^^ cm~'^ for all the giants except fi Vel, for which both HEG and MEG measurements 
suggest that Uc is compatible with the low-density limit {n^ ^ 10^^ cm~'^). 

One remarkable result among the giants is that the coronal density in the case of the 
supergiant Canopus appears to be high — of the order of 3 x 10^^ cm~^ indicated by Mg XI 
at temperatures ~ 10^ K. Consequently, the surface filling factor is tiny — a few millionths of 
the stellar surface. We discuss filling factors for the sample as a whole in more detail below. 
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4-2.3. Active Binaries 

Conspicuous among the active binaries is the clear indication of high densities — ric > 
10^^ cm~^ — from the Mg hne analysis for most systems. Most of the upper limits found cor- 
respond to HEG measurements that are generally affected by larger statistical uncertainties. 

Several previous analyses of spectra of the active stars of our sample exist, as summarized 
in Table 10. EUVE spectra of active binaries were analyzed by Sanz-Forcada et al. (2001, 
2002), who estimated plasma densities using Fe XIX- XXII lines (T~10 MK), and invariably 
found Tie ~ 10^^-10^^ cm~^. Other earlier EUVE results are summarized by Drake (2001) 
and Bowyer, Drake & Vennes (2000); Ness et al. (2002a, 2002b) presented an analysis of 
X-ray spectra obtained with the Chandra Low Energy Transmission Grating Spectrometer 
(LETGS) that provides a slightly higher sensitivity to the O VII lines than the HETG-FACIS- 
S combination. 

Some inconsistencies exist between our and earlier results for high temperature plasma 
(~ 10^ K). For example, Brickhouse & Dupree (1998) estimated > 10^^ cm~^ up to 
~ 10^^ cm~^ diagnosed from Fe XIX-XXII (T ~ lOMK) lines in EUVE spectra of 44 Boo. 
Again, our Si XIII, and possibly Mg XI, densities seem incompatible with the EUV estimates 
based on Fe lines, suggesting that the latter might be spuriously high. As previously men- 
tioned, the effect of lines blending with the density- sensitive diagnostics (e.g. Drake 2001), 
or an inaccurate theoretical Si XIII R ratio, might be responsible for this disagreement. 

4.2.4. General Trends 

Taking advantage of the large star sample analyzed, we investigated the presence of 
any specific trend related to the physical characteristics of the observed stars. We mainly 
concentrated on the results obtained from the Mg XI measurements, since this is the most 
complete sample. 

We first investigated trends in R ratio with stellar fundamental parameters such as effec- 
tive temperature and surface gravity. A correlation with surface gravity might be expected, 
for example, because surface gravity dictates the density and pressure that would charac- 
terise the plasma under hydrostatic equilibrium. The underlying stellar effective temperature 
might be expected to have a more subtle influence on coronal gas properties. In fact, we 
did not find any evident trend of the measured R ratios or derived densities with surface 
gravity or effective temperature. The most graphic illustration of this is when contrasting 
the giants, including the supergiant Canopus, and the dwarfs. R ratios are indistinguishable 
between these groups. 



- 19 - 



Perhaps more surprising, we also do not see obvious trends of the measured R ratio 
with rotation period or Rossby number. Instead, the only obvious trend in the R ratios is 
with stellar X-ray luminosity, and in the ratio of X-ray luminosity to bolometric luminosity, 
Lx/Lboi- The former trend is shown in Figure 13, where we show the O VII, Mg XI and 
Si XIII R ratios as a function of both Lx and the surface flux Fx. While there is arguably 
no discernible trend in any of the ratios with Fx, the Mg XI ratio shows a decrease in R 
ratio for Lx > 10^*^ erg s~^ or so. It is tempting to see the same trend in the Si XIII R 
ratio for Lx > 10'^^ erg s~^, but there are simply too few observations at the highest X-ray 
luminosities. The O VII R ratio is essentially the same for the whole range of surface fluxes 
and X-ray luminosities. 

The trend of declining Mg XI R ratio in more active stars becomes more clear when 
shown as a function of Lx/Lboi, as shown in Figure 16. We consider the derived densities 
for O VII and Mg XI as a function of Lx and Fx in Figure 15: a significant fraction of 
the measurements here formally yield upper limits to the plasma density and these tend to 
obscure the results that are more obvious in Figure 13 — it is only when the measured R 
ratio sample is taken together than the trends clearly emerge. Nevertheless, the trend of 
increasing density derived from Mg XI is still apparent at Lx > 10^^ erg s~^ or so; below 
this value the results are characterised more by upper limits. 

We emphasise that this trend of Mg XI R ratio and derived density with Lx is not 
determined by a possible selection effect related to the signal-to-noise ratio of the HETG 
spectra; in fact the relative error on the Mg XI R measurement, 5(R)/R, has no specific 
trend with Lx. 

4-2.5. Density correlated with Ne/Mg? 

In §3.3, we pointed out an apparent trend of measured Mg XI R ratio with Ne X/Mg XI 
line strength seen in Figure 9. We do not believe this trend is a result of residual blending; 
instead we suggest that we are seeing a trend of coronal abundance variation in which the 
Ne/Mg ratio is higher in stars characterised by higher Mg XI R ratios — higher density at 
temperatures of ~ 10^ K. Such a result might be expected based on the emerging trend 
of coronal abundance anomalies in active stars. The coronae of RS CVn-type binaries, in 
particular, exhibit high Ne/Fe abundances compared to their expected photospheric com- 
positions (see e.g. reviews by Drake et al. 2003, Audard et al. 2003). Since we see in this 
study good evidence for a trend of decreasing Mg XI R ratio with increasing Lx/i^boi, the 
correlation of Ne X/Mg XI line strength with Mg XI R ratio would be expected if the Ne/Mg 
abundance ratio increased with Lx/i^boi- While the Ne/Mg abundance ratio remains to be 
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investigated thoroughly in the stars of our sample, the evidence from analyses to date sum- 
marised by Drake (2003) and Audard (2003) does suggest that the relative Ne enhancement 
increases with stellar activity. 

4.3. Filling factors 

One of the most important aspects of density measurements is the associated insight 
into the structuring of the plasma and into the emitting volume. In the stellar application, 
we can also use the density information to investigate coronal filling factors. For coronal 
line emission we can write I^^^ = Jyn'^Gk{T,nc)dV, for the intensity of a line k emitted by 
optically thin plasma. In the isothermal approximation we have I^^^ = Gk(T,ne)EM where 
EM = nlV is the volume emission measure and Gk{T,ne) is the contribution function of 
the emission line k. If we further assume that the plasma is at the temperature of maximum 
formation of the line k, we end up with a lower limit for the EM of the plasma emitting that 
line. In turn, from EM we obtain an estimate of the surface coronal filling factor defined as 

V 1 _ EM 1 _ I^^' 1 1 
^ ~C'A~~^'A~ nlGk{T,n,) ' 2' 

where C is the scale height of the emitting plasma and Ai, is the surface area of the star. 
A possible choice of C is given by the hypothesis that the emitting coronal plasma is con- 
fined in hydrostatic loops smaller than, or comparable to, the pressure scale height. Under 
this assumption we can use the scaling laws of an hydrostatic loop model, such as that of 
Rosner, Tucker & Vaiana (1978; RTV) for which C ~ [T/(1.4 x 10^)]^/p. In this way, from 
the measured densities and assuming a temperature equal to the temperature of maximum 
formation of the line we can derive an estimate of the filling factors. 

In Fig. 17 we illustrate the surface filling factors, /, derived from Mg XI as a function 
of the rotation period and of the Rossby number. There is a sharp increase in filling factor 
at rotation periods of 3-4 days, below which the filling factor appears to saturate at values 
of several percent of the stellar surface, and above which there is an approximately linear 
trend between logarithmic rotation period and filling factor. A rotation period of 3-4 days 
also marks the general location of coronal "saturation" , and the saturation of filling factor 
here is not a surprise. We note that Canopus, whose rotation period is not accurately known 
(e.g. Decin et al. 2003), is missing from Fig. 17. Nevertheless, if its rotation period is of 
order several hundred days it would not be out of place in the relation between filling factor 
and period for P^ot > 4d. 

Figure 18 illustrates the filling factors derived from both Mg XI and O VII lines vs. 
X-ray luminosity and surface flux. In order to have a more complete sample and investigate 
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the relation with X-ray emission, we assumed an electron density of 2 x 10^°cm~^ for the 
stars whose O VII lines were not measurable. This is not unreasonable given that almost all 
the spectra of our active stars (with rare exceptions like the case of TW Hya) yield Uc close 
to this value; this estimate is also supported by analyses of Chandra LETGS and XMM- 
Newton observations of some of the sources of our sample for which we did not obtain O VII 
measurements, like 31 Com (Scelsi et al. 2004), Algol (Ness et al. 2002b), UX Ari (Ness et 
al. 2002a) and Prox Cen (Giidel et al. 2002). Nevertheless, it is notable that we find filling 
factors > 1 — a nonsensical value — for three of the active binaries (ER Vul, TY Pyx, and 
AR Lac) for which we assume the common density value, indicating that the true density 
at a few 10^ K in these stars must be higher than 2 x 10^°cm~^, or that the assumed scale 
height is too low. 

With some degree of scatter, the O VII filling factor, /ovii, is directly proportional to 
the X-ray surface flux. It is interesting that this filling factor begins to saturate — i.e. reach 
values between 0.1 and 1 — at a mean surface flux of ~ 10^ erg cm~^ s"^. While our data 
are sparse here, this saturation of /ovii seems to occur at surface flux values approximately 
one order of magnitude lower than the surface fluxes of the most active stars in the sample. 
The same value of surface flux, 10'' erg s~^ cm~^, was found for the areas of the solar surface 
covered with active regions by Withbroe & Noyes (1977). This suggests then, that at plasma 
temperatures of up to a few 10® K these "O VII saturated" stars are essentially covered in 
active regions possibly similar in nature to those of the Sun. 

The behaviour of the Mg XI filling factor is different. Firstly, there is a clear break in 
the fugXi-Fx relation at the same surface flux level, Fx ~ 10'^ erg s~^ cm~'^, as we see /ovii 
saturate. At lower Fx levels, the fugXi-Fx relation appears to have a shallower slope than 
that for /ovii vs. Fx though the level of scatter precludes a definitive statement. 

Finally, we examine the relative behavior of the two filling factors, /ovii and fugxi in 
Fig. 19. For /ovii ^ 0.2, the two factors are approximately proportional to one another. For 
/ovii ^ 0.2, fugXi increases much more sharply than /ovii. At higher activity levels (i.e. 
higher /) the coronae are preferentially filling with hotter plasma with respect to the cooler 
plasma emitting the O VII lines. 

4.4. Summary 

The density and filling factor results suggest a scenario that has been outlined before 
based on the emission measure distributions as a function of temperature of active dwarfs 
as compared to those of solar active regions by Drake et al. (2000; see also Giidel 1997): 
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going toward activity levels higher than the Sun the trend appears to be one of increased 
filling factor of solar-like active regions characterised by essentially the same plasma density 
of rie ~ 2 X 10^*^ cm~'^. For these active stars, there is an attendant hot, high density plasma 
with ~ 10^^ cm~^ whose contribution to the total X-ray increases toward higher activity 
stars. 

As the surface becomes "saturated" , or full, of active regions, with filling factors between 
0.1 and 1, the nature of the corona begins to change. This surface saturation occurs at X- 
ray luminosities well below the highest observed for the most active stars — by an order of 
magnitude or two. Drake et al. (2000) and Giidel (1997) suggested that interaction between 
adjacent active regions could give rise to increased fiaring behaviour which is manifest as 
the growth of the hotter 10^ K corona characterised in the Mg XI R ratio by higher plasma 
densities, and with a concomitant increase in the filling factor of this plasma. 

The quite different gas densities found for "high" and "low" coronal temperatures rein- 
forces the fact that the dominant emission at these temperatures cannot originate from the 
same structures: models seeking to explain active coronae need to account for the increase 
in gas pressure with increasing temperature. 

5. Conclusions 

We have investigated the density of coronal plasma at different temperatures through the 
analysis of the He-like triplet of O, Mg, and Si in high spectral resolution X-ray observations 
obtained with Chandra of a wide sample of active stars. As in previous works in the literature 
(e.g. Brickhouse 2002, Ness et al. 2003) the unprecedented Chandra spectrometer resolution 
proves to be vital for an accurate analysis: in this work we have especially emphasised the 
effect of lines blending with the Mg XI triplet lines and the need to take these into account 
in order to obtain reliable measurements of the R ratio. Ignoring these blends results in R 
ratios significantly different from values obtained when blends are accounted for. 

Our measured R and G ratios provide some insights into the propriety of recent theo- 
retical calculations of the relative line strengths from He-like ions. In the case of Si XIII, 
we find some evidence suggesting that the low density limit f /i ratios predicted by both 
APED and Porquet et al. (2001) are too low by 20-40% or so based on our observed values. 
Temperatures inferred from theoretical G ratios also tend to be lower than expected. We 
tentatively identify the cascade contributions from dielectronic recombination as the likely 
culprit for modelling inaccuracies. 

For most of the observed sources in which O VII yielded a useful R ratio measurement. 
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we find for plasma at temperatures of a few 10^ K densities of ~ 2 x 10^° cm~^. Mg XI lines 
formed closer to 10^ K yielded densities rie ~ 10^^ cm~'^ confirming the scenario outlined 
in the recent literature in which hotter (~ 6 — 10 MK) plasma is characterized by densities 
two orders of magnitude higher than the cooler plasma, and in which the latter has both 
temperature and density similar to those characterising solar coronal active regions. No 
star shows good evidence for Si XIII R ratios significantly below the theoretical low density 
limits; while this might partly be a problem associated with the theoretical ratios, which 
appear too low, it seems clear that there are no stars with electron densities significantly 
above 10^'^ cm~^. This result casts some doubt on earlier studies based on Fe XX- XXII 
lines seen in EUVE spectra which suggested the presence of densities exceeding 10^'^ cm~^ 
in active stellar coronae. 

We have found a distinct correlation between the measured Mg XI R ratio and X- 
ray luminosity, with lower R ratios, and by inference, higher densities, found in more X-ray 
luminous coronae. This trend is also visible as a function of the X-ray "production efficiency" , 
Lx/l-iboi- Instead, no correlations were found between R ratio and stellar parameters surface 
gravity and effective temperature, or with Rossby number. 

The surface coronal filling factors obtained from our density estimates are significantly 
smaller for the hotter plasma with T ~ 10^ K, with fugxi ~ 10~^ — 10~^, than for the 
cooler ~ 2-3 x 10^ K plasma characterized by /ovii values in the range 10~^ — 1. The 
remarkably small filling factors for hot plasma suggest the potential for studying X-ray 
rotational modulation using the hotter lines available in high resolution X-ray spectra. In 
this regard, we note that Marino et al. (2003) have recently found for the first time X-ray 
rotational modulation in a supersaturated star; Orlando et al. (2004) also show that such 
rotational modulation may be present in the hardness ratio. 

The O VII filling factor, /ovii, is found to be directly proportional to the X-ray surface 
fiux, and reaches values between 0.1 and 1 at a mean surface fiux of ~ 10^ erg cm~^ s~^ — the 
same surface fiux that characterises solar active regions. Instead, at higher temperatures the 
fugXi-Fx relation shows a break at the same surface flux level, Fx ~ 10'^ erg s~^ cm~^, after 
which it increases more strongly with increasing surface fiux. 

The density and filling factor results suggest the following picture: at plasma tempera- 
tures of up to a few 10^ K, stars with increasingly higher activity level become more and more 
covered with solar-like active regions with no ~ 2 x 10^'^ cm~^. At the same time, there is a 
similar growth in hotter plasma with temperatures close to 10'^ K, and with rie ~ 10^^ cm~^. 
This hotter plasma is likely associated with continuous superimposed fiaring activity associ- 
ated with the active regions, analogous to the solar case. For stars at the "O VII saturated" 
level, when X-ray surface fiuxes are similar to that of active regions on the solar surface. 



-24- 



the whole star is essentially covered by active regions. Further growth in X-ray emission at 
higher activity levels still is in the form of increased flaring, probably arising because of the 
crowding of active magnetic structures on the stellar surface. This activity gives rise to a 
growth in the surface coverage of 10^ K plasma at higher densities of n^. ~ 10^^ cm~^. This 
picture is essentially that described earlier by Giidel (1997) and Drake et al. (2000). 
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Ministero dell'Istruzione, dell'Universita e della Ricerca and by Agenzia Spaziale Italiana. 
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Table 1. List of stellar parameters. 



Source HD Spcctr. 

Typo 



d " 
[pc] 



R^/R(r 



M^/Mq 



log(Lboi) Porb 
[erg/ sec] [days 



[days] 



AU Mic 


197481 


MIV 


9.9 


0.56 ^ 


0.59 1 


3730 27 




1.44 37 


32 


.48 


39 




4.854 "13 


0. 


024 


Prox Cen 




M5Vc 


1.3 


0.16 ^ 


0.11 =0 


2700 28 




1.97 


30 


.58 






83.5 " 






EV Lac 




M4.5V 


5.1 


0.41 ^ 


0.34 3 


3300 29 




1.57 


31, 


.62 


3 




4.376 *5 






AB Dor 


36705 


KOV 


15 


1. ^ 


0.76 " 


5000 39 




0.80 38 


33 


.18 


40 




0.51479 


0. 


,0038 


TW Hya 




K8Vc 


56 


1. = 


0.7 ^ 


4150 31 




0.7 38 










2.9 " 


0. 


,019 




223460 


Gl III 


135 


13.6 


2.9 21 


5110 21 




0.81 6 


35, 


.43 






23.25 6 


0. 


19 


31 Com 


111812 


GO III 


94 


8 ^ 


2.96 22 


5320 21 




0.68 


35, 


.41 


7 




> 5 "8 


0. 


,23 


/3 Ceti 


4128 


KOIII 


29.4 


15.1 ^ 


3.2 21 


4840 21 




1.2 ^ 


35 


.77 






> 255 21 


0. 


,27 


Canopus 


45348 


FOII 


95.9 


53 


13 23 


7350 23 




0.15 38 


37 


.70 


41 




> 298 23 






IJ. Vel 


93497 


G5III/... 


35.5 


13 9 


3 9 


4862 32 




0.90 


35, 


.62 






> 103 32 


0. 


,38 


Algol 


19356 


B8V/K1IV 


28 


2.9/3.5 ^° 


3.7/0.8 2-1 


11400/5300 


24 




34, 


.47 


42 


2.8 19 


2.8 19 


0. 


024 


ER Vul 


200391 


G0V/G5V 


50 


1.07/1.07 11 


1.10/1.05 


6000/5619 


33 


0.68 


33, 


.89 


0.6981 


0.6942 


0. 


,0077 


44 Boo 


133640 


G1V/G2V 


13 


0.87/0.66 1^ 


0.98/0.55 12 


5300/5035 


12 


0.65 


33, 


.82 


42 


0.267 




0. 


,0036 


TZ CrB 


146361 


GOV/GOV 


22 


1.1/1.1 " 


1.12/1.14 


6000/6000 


13 


0.51 


33 


.88 


1.1398 


1.1687 


0. 


,016 


UX Ari 


21242 


G5V/K0IV 


50 


1.11/5.78 l* 


> 0.63/ > 0.71 


.../4800 3-1 


0.91 


34 


.49 


42 


6.4379 


6.4379 


0. 


042 


5 UMa 


98230 


GOV/GOV 


7.7 


0.95/... 1^ 


0.92/... 25 


5650 2= 




0.59 


33 


.35 


42 


3.9805 




0. 


044 


II Pog 


224085 


K2V/... 


42 


3.4/... 


0.8/... 


4600 16 




1.01 


34, 


.17 


49 


6.7242 


6.718 


0. 


,056 


A And 


222107 


G8III/... 


26 


7.4/... 1^ 


0.65/... 26 


4800 26 




1.01 


34 


.80 




20.521 


53.952 


0. 


,31 


TY Pyx 


77137 


G5IV/G5IV 


56 


1.59/1.68 


1.22/1.20 1^ 


5400/5400 


IS 


0.72/0.76 


34 


.26 


42 


3.1986 


3.32 


0. 


,25 


AR Lac 


210334 


G2IV/K0IV 


42 


1.8/3.1 11 


> 1.30/ > 1.30 


.../5100 35 


0.72 


34, 


.49 


42 


1.9832 


1.9832 


0. 


,018 


HR 1099 


22468 


G5IV/K1IV 


29.0 


1.3/3.9 11 


1.1/1.4 


5400/4800 


36 


0.92 11 


34, 


.31 


42 


2.8377 


2.841 


0. 


,0076 


IM Peg 


216489 


K2III-II/... 


96.8 


13/... 19 


1.5/... 19 


4450 19 




1.12 


35 


.35 


42 


24.65 


24.39 


0. 


,036 
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Table 2. Parameters of the HETG observations. 



Source 


Obs ID 


Lx ^ 


Lx^ 


Fx ' 


^exp 






[erg/sec] 


[erg/sec] 


[10^erg/cm2/sec] 


[ksec] 


AU Mic 


17 


1.05-1029 


1.29-1029 


60.7 


58.8 


Prox Cen 


2388 


3.33-1026 


4.56-1026 


2.93 


42.4 


EV Lac 


1885 


2.59-1028 


3.19-1028 


31.1 


100 


AB Dor 


16 


6.83-1029 


8.26-1029 


136 


52.3 


TW Hya 


5 


1.15-10^" 


1.30-103" 


214 


47.7 


HD 223460 


1892 


6.31-10^1 


5.96-1031 


52.9 


95.7 


31 Com 


1891 


4.85-1030 


5.15-103" 


13.2 


130.2 


(3 Ceti 


974 


2.16-103" 


2.62-103" 


1.89 


86.1 


Canopus 


636 


2.48-103" 


2.94-103" 


0.17 


94.6 


/i Vel 


1890 


1.19-103" 


1.41-103" 


1.37 


19.7 


Vel 


3410 


1.17-103" 


1.45-103" 


1.4 


57 


Algol 


604 


8.50-103" 


8.88-103" 


119 


51.7 


ER Vul 


1887 


2.46-103" 


2.78-103" 


399 


112.0 


44 Boo 


14 


3.59-1029 


4.56-1029 


172 


59.1 


TZ CrB 


15 


2.99-103" 


3.34-103" 


453 


83.7 


UX Ari 


605 


5.69-103" 


7.26-103" 


35.7 


48.5 


C UMa 


1894 


8.78-1028 


1.23-1029 


22.4 


70.9 


II Peg 


1451 


1.56-1031 


1.76-1031 


250 


42.7 


A And 


609 


1.34-103" 


1.85-103" 


5.55 


81.9 


TY Pyx 


601 


4.71-103" 


5.10-103" 


573 


49.1 


AR Lac 


6 


5.21-103" 


5.60-103" 


284 


32.1 


AR Lac 


9 


5.61-103" 


6.30-103" 


319 


32.2 


HR 1099 


62538 


7.85-103" 


1.05-1031 


113 


94.7 


IM Peg 


2527 


2.75-1031 


2.79-1031 


27.1 


24.6 


IM Peg 


2528 


2.17-1031 


2.30-1031 


22.3 


24.8 


IM Peg 


2529 


1.86-1031 


1.97-1031 


19.2 


24.8 



'^relative to the HEG range : 1.5-15A 
•"relative to the MEG range: 2-24A 
'^from Lx obtained from MEG spectra 

Note. — In order to calculate the X-ray surface flux, Fx, for binary 
systems, we used the radius i?* of the component which probably is the 
stronger X-ray source, e.g. the KlIV component for Algol. 



Table 3. Theoretical wavelength of the He- like triplet lines. 



ion 




wavelength (A) 






resonance (r) 


intercomhination (i) 


forbidden (f) 




ls2 iSo - ls2p 


ls2 ISo - ls2p ^V2S 


ls2 iSo - ls2s ^Si 


Si XIII 


6.6479 


6.6850 / 6.6882 


6.7403 


Mg XI 


9.1687 


9.2282 / 9.2312 


9.3143 


VII 


21.6015 


21.8010 / 21.8036 


22.0977 



Table 4. Neon H-like Lyman series. 



transition 


Lya: 2->l 


Ly/3: 3^1 


4^1 


5^1 


6^1 


7^1 


8^1 


9^1 


10^1 


A (A) 


12.134 


10.239 


9.708 


9.481 


9.362 


9.291 


9.246 


9.215 


9.194 


r 


1.0000 


0.1900 


0.0697 


0.0335 


0.01874 


0.01157 


0.0076 


0.0053 


0.0038 



oscillator strength normalized to the /-value of the Lya transition. 
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Table 5. Mg XI triplet lines and blending components. 



^obs 

[A] 


Ion 




transition 


Emissivity^ 
[ph cm^s"^] 


9.1689 


Mg XI 


r 


ls2p -> ls2 




1. 10x10-16 


9.1900 


Fe XXI 




Is22s2p2 4p3/2 ^ Is22s22p2 Sp^ 


1.14x10-1^ 




Fe XX 








3.29x10-18 




Ne X 


10^1 




1x10-18 


9.2187 


Fe XX 




2s22p2(3p)5d 2p_,^^ ^ 


2s^2p^ 2^3/2 


1.63x10-18 




No X 


9^1 






1.4x10-18 


9.2304 


Mg XI 


i 


ls2p ^ ls2 15*0 




1.60x10"", 2.23x10-18 


9.2467 


Nc X 


8-^1 






2.1x10-18 


9.2882 


Ne X 


7^1 






3.1x10-18 




Fe XX 




2s22p2(3p)5rf 4p^/^^ 




2.82x10-18 




Fe XXII 




ls^2s2p{^P)Ad 4D3/2 - 


-> Is22s2p2 2^)3^2 


1.37x10-18 


9.3144 


Mg XI 


/ 


Ls2.s ^S*! ^ ls2 15-0 




5.31x10-1^ 



^From the fit to the HEG coadded spectrum of /? Ceti, TZ CrB, AR Lac, and HR 1099 shown in 
Fig. 7. 

I'Peak hne emissivity from the APED database for an electron density rie = 1 x 10i° cm~^ , 
and assuming the chemical composition of Anders & Grevcsse (1989). The excitation rate for a 
Maxwellian distribution of electron velocities for higher n series in hydrogenic ions is essentially 
proportional to the /-value of the transition (e.g., Lang 1999); since emissivities for the Ne Lyman 
scries lines do not appear to be available for upper level ri > 5, the emissivities listed here are the 
APED emissivity of the 5 — > 1 transition, scaled in proportion to the transition /-values oscillator 
strengths (see Tab. 4). 
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Table 6. Measured Line fluxes (in 10 ^ photons cm ^ sec ^) with la errors. 



Source grating flux 









S 


i XIII 












Mg XI 






O VII 




r 




i 






/ 




r 


i 


/ 


r 


i 


/ 


AU Mic 


HEG 


24 ± 5 




7 ± 


4 




19 ± 


4 
















MEG 


27.3 ± 2.8 


6, 


.1 ± 1. 


,9 


16, 


5 ± 2. 


1 


19.4 it 1.7 


5.0 it 1.8 


11.9 it 1.7 


180 it 35 


48 it 29 


150 ± 40 


Prox Ccn 


HEG 






























MEG 


7.5 ± 2.2 


2 


.0 ± 1. 


,8 


4 


2 ± 1. 


7 


5.2 ± 1.5 


2.2 ± 1.5 


6.5 it 1.5 








EV Lac 


HEG 


31 ± 4 


6, 


.6 ± 2. 


4 


20 


5 it 2. 


8 


18 it 4 


6 it 3 


9 ± 3 










MEG 


29.3 ± 2.0 


7 


.8 ± 1. 


,4 


17 


,3 it 1. 


5 


18.9 it 0.6 


6.3 it 0.7 


9.3 it 0.6 


207 it 26 


78 it 20 


117 ± 25 


AB Dor 


HEG 


56 ± 7 




13 ± 


5 




45 ± 




66 it 9 


15 it 7 


24 it 7 










MEG 


49 ± 4 


9, 


.8 ± 2. 


,4 




33 it 


3 


48.9 it 1.5 


13.8 it 1.5 


37.3 it 1.5 


276 it 45 


100 it 40 


176 ± 50 


TW Hya 


HEG 






























MEG 


2.8 ± 1.4 


1 


.8 ± 1. 


2 


3, 


,0 it 1. 


2 








90 it 30 


73 it 30 


< 4 


HD 223460 


HEG 


32 ± 4 




10 ± 


3 




24 ± 


3 
















MEG 


24.5 ± 2.2 


9, 


.6 ± 1. 


,8 


17, 


2 ± 1, 


,9 


19.4 ± 0.9 


5.6 ± 1.0 


10.7 ± 1.0 


36 ± 17 






31 Com 


HEG 


19.0 ± 2.6 


4, 


.1 ± 1. 


6 


9 


3 it 1. 


8 


13 it 3 


4.0 it 2.2 


6.1 it 2.5 










MEG 


21.1 ± 1.5 


3, 


.6 ± 0. 


,9 


12 


,6 it 1. 


1 


15.2 it 0.5 


3.5 it 0.5 


7.1 ± 0.6 


28 it 13 






P Get 


HEG 


81 ± 6 




21 ± 


3 




66 it 


5 


113 ± 7 


26 ± 5 


78 ± 7 










MEG 


81 ± 3 


20, 


,6 ± 2. 


1 


49, 


,9 it 2. 


5 


108 it 2.0 


28.7 it 1.2 


64.5 it 1.2 


100 it 24 






Canopus 


HEG 


8.6 ± 1.7 


1 


.9 ± 1. 


7 


6 


2 it 1. 


3 
















MEG 


8.2 ± 1.2 


1 


.5 ± 0. 


,8 


6, 


7 it 1. 


,0 


9.0 it 0.7 


2.6 it 0.7 


5.0 it 0.7 


31 it 16 






fl Vel 


HEG 


33 ± 5 


10 ± 3. 


5 




20 it 


4 


57 ± 7 


10 ± 5 


31 ± 6 










MEG 


30.4 ± 2.4 


9, 


.9 ± 1. 


7 


19, 


,3 it 1. 


,9 


45.0 it 1.1 


7.6 ± 1.16 


23.9 ± 1.2 


45 it 26 






Algol 


HEG 


133 it 10 




36 ± 


7 




85 it 


8 


126 ± 13 


31 ± 10 


78 ± 12 










MEG 


125 ± 7 




31 ± 


5 




76 it 


4 


102 it 2.8 


24.0 it 2.3 


57.6 it 2. 1 


200 it 33 






ER Vul 


HEG 


22 ± 3 


6, 


.7 ± 2. 


,3 


13, 


,8 it 2. 


5 


33 ± 4 


8 ± 3 


19 ± 4 










MEG 


24.1 ± 1.7 


4, 


.6 ± 1. 


1 


15, 


,0 it 1. 


,3 


22.9 it 0.6 


5.59 it 0.6 


15.2 it 0.6 


60 it 19 






44 Boo 


HEG 


50 ± 7 




11 ± 


4 




38 it 


5 


95 ± 10 


21 ± 7 


52 ± 9 










MEG 


44 ± 3 


8, 


.9 ± 2. 


,4 


28, 


,9 it 2. 


,6 


63.1 it 2.5 


14.8 it 1.2 


38.6 it 1.2 


280 it 50 


100 it 40 


170 ± 50 


TZ CrB 


HEG 


123 ± 7 




20 ± 


4 




72 it 


5 


176 ± 9 


45 ± 6 


89 it 8 










MEG 


129 ± 4 


25, 


,0 ± 2. 


7 




73 it 


4 


152 ± 4 


33.0 it 2.8 


78 it 2.0 


280 it 40 


99 it 27 


230 ± 40 


UX Ari 


HEG 


33 ± 7 




14 ± 


5 




24 it 


5 


31 ± 7 


15 ± 7 


22 ± 7 










MEG 


31 ± 3 


7, 


.3 ± 2. 


,4 


24, 


,0 it 2. 


7 


27.5 it 1.8 


12.4 it 1.3 


17.9 ± 1.3 


150 it 40 


27 it 27 


60 ± 30 


5 UMa 


HEG 


41 ± 5 




13 ± 


4 




27 it 


4 


76 ± 7 


19 ± 5 


55 ± 7 










MEG 


46 ± 3 


8, 


.4 ± 1. 


,8 


24, 


7 it 2. 


1 


68.1 it 2.1 


18.7 it 0.9 


42.4 it 1.2 


360 it 40 


84 it 29 


200 ± 40 


II Peg 


HEG 


95 ± 10 




27 ± 


7 




50 it 


7 


68 ± 11 


32 ± 10 


46 ± 11 










MEG 
















60.3 it 2.2 


20.3 it 1.8 


34.3 it 2.1 


255 it 60 


127 it 50 


210 ± 60 


A And 


HEG 


61 ± 5 




14 ± 


3 




41 it 


4 


93 ± 7 


15 ± 4 


45 ± 5 










MEG 


58 ± 3 


7, 


.6 ± 1. 


7 


31 


6 it 2. 


2 


84.5 ± 0.8 


16.8 it 0.7 


40.3 ± 0.8 


135 it 27 


36 it 22 


88 ± 28 


TY Pyx 


HEG 


29 ± 6 




8 ± 


4 




22 it 


5 
















MEG 


28 ± 3 


7, 


.7 ± 2. 


2 


19, 


4 it 2. 


5 


29.1 it 1.4 


6.7 it 1.5 


15.8 it 1.4 


77 it 30 






AR, Lac 


HEG 


61 ± 6 




19 ± 


5 




47 it 


5 


64 ± 9 


23 ± 7 


44 ± 9 










MEG 


51 ± 3 


10, 


,9 ± 2. 


,3 


30, 


,9 it 2. 


,6 


69.0 it 1.6 


14.8 it 1.4 


36.8 ± 1.5 


160 it 40 






HR 1099 


HEG 


83 ± 6. 




23 ± 


4 




64 it 


5 


108 ± 8 


26 ± 6 


55 ± 7 










MEG 


102 ± 4 


15, 


7 ± 2. 


,6 




59 it 


3 


97 ± 4 


23 ± 2 


52 ± 2 


400 it 40 


83 it 27 


240 ± 40 


IM Peg 


HEG 


34 ± 6 




8 ± 


6 




22 it 


5 
















MEG 


37 ± 3 


9, 


.3 ± 2. 


5 


17, 


,6 it 2. 


4 


25.7 it 1.6 


7.1 it 1.7 


13.5 ± 1.7 


30 it 16 







^r, resonance line, i, intercombination lines, f, forbidden line, corresponding to the transitions described in Table 3. 
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Table 7. Line flux measurements. 



Source flux 



Mg XII Ne X O VIII Fe XVIII Fe XXI 





Lya 
8.419A 


Lya 
12.132A 


Lya 
18.967A 


14.210A 


12.292A 


AU Mic 


22 ±5 


324 ± 13 


790 ± 40 


59± 11 


30 ±6 


Prox Ccn 




60 ±7 


220 ± 30 






hj\ Lac 


15.0 ± 1.8 


230 ± 8 


f^'TO 1 on 

d73 ± 29 


61 ± 7 


23 ± 3 


AB Dor 


64 ±3 


683 ± 19 


1470 ± 50 


201 ± 15 


169 ± 12 


TW Hya 




85 ±8 


240 ± 30 






HD 223460 


56 ±3 


129 ±6 


168 ± 18 


54 ±8 


50 ±5 


31 Com 


22.3 ± 1.7 


59 ±4 


97± 12 


37 ±5 


50 ±4 


P Get 


100 ±3 


397 ± 11 


596 ± 29 


403 ± 14 


136 ± 7 


Canopus 


8.5 ± 1.0 


55 ±4 


109 ± 15 


33 ±5 


15.6 ±2.8 


^l Vol 


39 ±3 


143 ± 78 


263 ± 26 


162 ± 13 


49 ±5 


Algol 


207 ±4 


872 ± 23 


1180 ±50 


259 ± 18 


258 ± 12 


ER Vul 


36 ±3 


188 ± 7 


262 ± 20 


84 ±7 


53 ±4 


44 Boo 


60 ±4 


541 ± 14 


1310 ±50 


158 ± 15 


75 ±8 


TZ CrB 


188 ±3 


1059 ± 16 


1830 ±50 


440 ± 19 


269 ± 10 


UX Ari 


48.1 ±2.4 


743 ± 20 


910 ±50 


76± 12 


58 ±7 


^ UMa 


50 ±3 


371 ±11 


1160 ±40 


173 ± 12 


44 ±5 


II Peg 


100 ±4 


1181 ±28 


1930 ±60 


106 ± 15 


105 ± 11 


A And 


116 ±3 


553 ± 12 


950 ± 30 


119± 10 


95 ±6 


TY Pyx 


60 ±4 


294 ± 13 


360 ± 30 


78± 12 


71 ±8 


AR Lac 


89 ±4 


632 ± 16 


820 ± 40 


161 ± 15 


126 ± 10 


HR 1099 


137 ±4 


1705 ±21 


2860 ± 60 


198 ± 12 


207 ±9 


IM Peg 


46 ±4 


334 ± 11 


340 ± 30 


76± 12 


61 ± 7 



fluxes, in 10 ^ photons cm ^ sec ^, from MEG spectra, with Icr errors. 
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Table 8. Densities. 



Source^ rio (cm 

Mg XI O VII 

HEG MEG MEG 



AU Mic 


- 




< 5.62 


1012 


< 5.62 


10" 


Prox Cen 


_ 




< 5.6 • 


1012 


_ 




EV Lac 


5.6if 


1012 


^•"-0.6 


•1012 


5.612-5 


•lOlo 


AB Dor 


< 3.2 • 


10" 


< 1.8 • 


1012 


5.6tl.4 


•lOlo 


TW Hya 










> 7.8 • 


10" 


HD 223460 






Q 9+2.5 


•1012 






31 Com 


5.6lf 


10^2 


Q 9+2.5 


•1012 






/3 Get 


< 1.8- 


1012 


1 C+0.6 


•1012 






Canopus 






"^•2-2.2 


•1012 






At Vel 


< 5.6 • 


10^2 


< 1.0 • 


1012 






Algol 


< 5.6 • 


1012 




•1012 






ER Vul 


< 5.6 • 


1012 


< 1.8 • 


1012 






44 Boo 


< 5.6 • 


1012 


1 0+°-^ 

^•^-0.8 


•1012 


1 0+^ 


•IQi" 


TZ CrB 


q 9+2.5 


•1012 


1 C+0.6 


•1012 




•lOio 


UX Ari 


5.6^f 


1012 




1012 






CUMa 


1.8+30 


10" 


q 9+0.6 
0.^_1.4 


•1012 




•lOio 


II Peg 




•1012 




•1012 


3.2^1, 


•lOio 


A And 


< 3.2 • 


1012 




•1012 


< 1.0 • 


10" 


TY Pyx 








•1012 






AR Lac 


3-2^2.2 


•1012 


< 1.8 • 


1012 






HR 1099 


q 9+2.5 
■^•^-2.2 


•1012 


1 c+0.6 


•1012 


1 0+2-2 


•lOi" 


IM Peg 






q 9+2.5 
0.^„1 4 


•1012 







'^Horizontal lines divide the different classes of source into 
single d-warfs (top), single giants (middle), and active binaries 
(bottom). 
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Table 9. Comparison of derived densities with previous work: single stars. 



This work Previous works 



Source^ 


logT^ 
[MK] 


log Tic 


lon'^ 


[MK] 


log 71c 

[cm~^] 


lon'^ 


data'^ 


Ref. 


AU Mic 


6.8 
6.3 


< 12.75 

< 11.75 


MgXI 
OVII 


7 

6.3 
6.3 
5.3 


- 13 

10.0 - 10.3 (flaring) 
< 9 (quiescent) 
10.7-11 


FeXXI 

OVII 

OVII 

OV, OIV, SIV 


D 
C 
C 

D,E,F 


1 

2 
2 

3,4 


Prox Cen 


6.8 


< 11.75 


MgXI 


6.3 
6.3 


~ 11.6 (flaring) 
~ 10.3 (quiescent) 


OVII 
OVII 


C 

c 


5 

5 


AB Dor 


6.8 
6.3 


< 12.25 


MgXI 
OVII 


7 

6.3 


12 - 13 
~ 10.8 


FeXXI-XXII 
OVII 


A,C,D 
A,C 


6,7 
8 


TW Hya 


6.3 


> 11.9 


OVII 


6.3 


> 12 


OVII 


A,C 


9,10 


HD 223460 


6.8 


12.5 ±0.25 


MgXI 


6.3-7 




OVII, NelX, SiXII 


C 


11 


31 Com 


6.8 


12.5 ±0.25 


MgXI 


6.3 

5.2 


9.8 - 11.3 
10.0-10.6 


OVII 
OIV 


C 

F 


12 
13 


P Get 


6.8 


12.25^Q 2 


MgXI 


7 

5.2 


11.9-13.1 
9.7-10.3 


FeXIX-XXII 
OIV 


D 
F 


6 

13 


fi Vel 


6.8 


< 12.25 


MgXI 








D 


14 



''Horizontal lines divide the different classes of source into single dwarfs (top) and single giants (bottom). 

''Temperature of maximum formation of the ion emitting the density sensitive lines used for the density diagnos- 
tics. 

'^lon emitting the density sensitive lines used for the diagnostics. 

^A: C/iandra-HETGS; B: Chandra-LETGS; C: XMM-Newton-RGS: D: EUVE; E: FUSE; F: HST. 



References. — (1) Monsignori Fossi et al. (1996); (2) Magee et al. (2003); (3) DelZanna, Landini & Mason (2002); 
(4) Redfield et al. (2002); (5) Giidel et al. (2002); (6) Sanz-Forcada et al. (2002); (7) Sanz-Forcada et al. (2003a); 
(8) Sanz-Forcada et al. (2003b); (9) Kastner et al. (2002); (10) Stelzer & Schmitt (2004); (11) Gondoin (2003); (12) 
Scelsi ct al. (2004); (13) Ayres et al. (1998); (14) Ayres et al. (1999). 
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Table 10. Comparison of derived densities with previous work: active binaries. 



This work Previous works 



Source logT^ logric Ion'' logT^ log rie Ion'' data'^ Ref. 

[MK] [cm-3] [MK] [cm-^] 



Algol 
ER Vul 

44 Boo 

TZ CrB 
UX Ari 

C UMa 
II Peg 



6.8 

6.8 

6.8 
6.3 

6.8 
6.3 

6.8 



6.8 
6.3 

6.8 
6.3 



+0.1 
0.5 



A And 6.8 
AR Lac 6.8 

HR 1099 6.8 
6.3 



12.3 
< 12.3 

1 0+0.3 
-^^-0.7 
1Q+0.8 
^^-0.2 



12.3 



+0.1 
-0.1 



13 



-0.03 
-0.3 



+0.1 



1 9 ^ 

10.31°:^ 
10.51;;^ 



12 3+°-^ 
12.5 ±0.5 



12 3+°-^ 

lo.otU:? 



MgXI 


6.1-6.3 


- 10.5 


NVI, OVII 


B 


1,2 


MgXI 


7 


12.2-12.8 


FeXIX-XXI 


D 


3 


MgXI 


7 


12- 14 


FeXIX-XXI 


D 


4 


OVII 












MgXI 


6.6-7 


12.0-12.5 


NelX, MgXI, SiXII 


A 


5 


OVII 


6.3 


10.0-10.8 


OVII 


A 


5 


MgXI 


7 


12.1 - 12.9 


FeXIX-XXI 


D 


6 




6.6 


11.1 - 11.7 


NelX 


B 


1 




6.1 


9.8-10.8 


NVI 


B 


1 


MgXI 


7 


11.8-12.3 


FeXIX-XXI 


D 


3 


OVII 












MgXI 


7 


12.4-13.4 


FeXIX-XXI 


D 


6 


OVII 


6.3 


10.6-11.6 


OVII 


A 


7 




6.8 


12.8-13.8 


MgXI 


A 


7 




6.6 


11 - 12 


NcIX 


A 


7 


MgXI 


7 


11.7-12.4 


FeXIX-XXI 


D 


6 


MgXI 


6.8 


^ 12.8 


MgXI 


A 


8 




6.3 


9-12 


OVII 


A 


8 


MgXI 


6.8 


11 - 12 


MgXI 


B 


1 


OVII 


6.6 


< 11 


NelX 


B 


1 




6.3 


10.1 - 10.6 


OVII 


A,B 


9,1 




6.3-7 


< 11 


MgXI, SiXIII 


A 


9 




7 


< 12 - 13 


FeXXI 


D 


9 



'^Temperature of maximum formation of the ion emitting the density sensitive lines used for the 
density diagnostics. 

''Ion emitting the density sensitive lines used for the diagnostics. 

'^A: Chandra-UETGS] B: C/iandra-LETGS; C: XMM-Newton-RGS; D: EUVE. 

References. — (1) Ness et al. (2002a); (2) Ness et al. (2002b); (3) Sanz-Forcada et al. (2003a); (4) 
Brickhouse & Dupree (1998); (5) Osten et al. (2003); (6) Sanz-Forcada et al. (2002); (7) Huenemoerder 
et al. (2001); (8) Huenemoerder et al. (2003); (9) Ayres et al. (2001). 
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Fig. 1. — HEG (black) and MEG (grey) spectra of the single stars of the sample. 
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Fig. 2. — HEG (black) and MEG (grey) spectra of the analyzed active binary systems. 



- 35 - 




Fig. 3. — Lightcurves, obtained as the sum of total counts of the HEG and MEG dispersed 
spectra over temporal bins of 100 sec, for the single stars of the sample. In each plot the 
initial time in MJD is indicated. 
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Fig. 4. — Lightcurves, as in Figure 3, for the RS CVn systems. 




Fig. 5. — P Ceti MEG spectrum. The analyzed hues are marked. 
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Fig. 6. — Theoretical R values from APED (Smith et al. 2001), at the temperature of 
maximum formation of each triplet, plotted vs. Uc for the three analyzed He-like triplets. 
The gray curves correspond to the theoretical values from Porquet et al. (2001). 
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Fig. 7. — Mg XI He-like triplet spectral region: Top- HEG spectrum of HR 1099: the best 
fit model, with only the continuum and the He-like triplet components, is superimposed to 
the data points with their associated error bars. Middle- coadded HEG spectra of four of 
the sources with the highest S/N: /3 Ceti, TZ CrB, AR Lac and HR 1099. The best fitting 
model with the blending components is superimposed to the spectrum. Bottom- Same as in 
the middle panel but for MEG spectra. 
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Fig. 8. — Values of R {upper panels) and G [lower panels) ratios for Mg XI, obtained taking 
into account the blending with the Fe and Ne lines, plotted vs. the values derived not 
considering the blending. The dashed lines mark the locus of equality. Different symbols are 
used for different source types: circles for single dwarfs, diamonds for giants, triangles for 
binary systems. 
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^multiple systems 



012345 012345 1 10 

FeXVIIl(l4.21A)/MgXI(r) FeXXl(l 2.29A)/MgXI(r) NeX(Lya)/MgXl(r) 

Fig. 9. — Investigation of possible correlations of the R ratio with Fe or Ne intensity in order 
to check for blending effects on Mg triplet lines. The dashed line marks the limiting value 
of the R ratio; lower values of R correspond to higher density (see Fig. 6). Data points 
from HEG measurements are shifted by +10% along the x-axis, to distinguish the error bars 
corresponding to HEG and MEG measurements for the same source, otherwise overlapping. 
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Fig. 10. — The Si XIII He-like triplet, as observed with the HEG, and examples of spectra 
and model fits for the three different type of sources: single stars {left), giants (center) and 
binary systems (right). Data points are shown with their associated error bars and a solid line 
superimposed to the data correspond to the best fitting model; we also show the components 
of the model, i.e. the continuum emission and the triplet emission lines. 
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Fig. 11. — Examples of analyzed spectral regions, as in Figure 10, for MEG spectra in the 
Mg XI triplet region, for single stars (left), giants (center) and binary systems (right). 
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Fig. 12. — Examples of analyzed spectral regions, as in Figure 10, for MEG spectra of the 
O VII triplet region, for single stars (left), and binary systems (right). None of the spectra 
of giant stars yielded significant measurements. For better readability we omitted the error 
bars in the TW Hya spectrum. 
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Fig. 13. — Measured R ratios for the He-like emission triplets, Si XIII {left), Mg XI {center) 
and O VII {right), from both HEG (filled) and MEG (empty) data, vs. the stellar X-ray 
surface flux {upper panels) and the X-ray luminosity {lower panels) for all the analyzed 
stars. The dashed line marks the limiting value of the R ratio; lower values of R correspond 
to higher density (see Fig. 6). Symbols as in Figure 9. 
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Fig. 14. — Comparison between the temperatures derived from the G ratio of He-hke triplet 
hnes, and the temperatures derived from the ratio of the Lya line of the H-hke ion and the 
resonance line of the He-like ion, for Mg [left) and O [right). The dashed hnes mark the 
locus of equal values. Symbols as in Figure 9. 
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Fig. 15. — Plasma densities derived from the R ratio of He-like triplets of Mg {upper panels) 
and O [lower panels) as a function of the surface X-ray flux [left) and vs. the X-ray luminosity 
(right). Error bars compatible with the lower or the upper limit of the density range of 
sensitivity of the R ratio are represented as down and up arrows, respectively. The upper 
limit symbols mark densities compatible with the low-density limit. Symbols as in Figure 9. 
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Fig. 16. — R ratios from the Mg XI triplet plotted vs. Lx/Lboi- The dashed line marks the 
limiting value of the R ratio; lower values of R correspond to higher density (see Fig. 6). 
Symbols as in Figure 9. 
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Fig. 17. — Surface filling factors estimated from the measured densities, vs. the rotation 
period (left) and vs. the Rossby number (right) of the sources. Symbols as in Figure 9. 
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Fig. 18. — Surface filling factors derived from Mg XI (top) and O VII {bottom) lines vs. 
X-ray surface flux (left) and the X-ray luminosity (right) of the sources. Symbols as in 
Figure 9. The gray symbols in the bottom panels mark the sources for which we assumed 
ne = 2 X 10^° cm-3. 
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Fig. 19. — Mg XI surface filling factors vs. O VII surface filling factors. Symbols as in 
Figure 9. The dashed lines superimposed to the data represent fugxi oc /ovii and fugxi oc 
/ovii- 
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